Creep motions of flux lines in type II superconductors with point-like defects 
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We simulated the creep motions of flux lines subject to randomly distributed point-like pinning 
centers. It is found that at low temperatures, the pinning barrier U defined in the Arrhenius-type 
v — F characteristics increases with decreasing force U{F) oc as predicted by previous theories. 

The exponent jj, is evaluated as 0.28 ± 0.02 for the vortex glass and fi ~ 0.5 ± 0.02 for the Bragg 
glass (BrG). The latter is in good agreement with the prediction by the scaling theory and the 
functional-renormalization-group theory on creep, while the former is a new estimate. Within BrG, 
we find that the pinning barrier is suppressed when temperature is lifted to approximately half of 
the melting temperature. Characterizations of this new transition at equilibrium are also presented, 
indicative of a phase transition associated with the replica-symmetry breaking. 
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PACS numbers: 74.25.Qt; 74.25.Dw 

Introduction - An elastic medium immersed into random 
environment is a generic system modeling a wealth of 
physical situations. In spite of the considerable efforts 
expended and impressive progress, many crucial issues 
of the underlying physics remain unresolved, and under- 
standing the statics and dynamics of such systems is one 
of the major challenges of condensed matter physics. 

Intensive activities have taken place in vortex states of 
type II superconductors since the discovery of high-Tc su- 
perconductivity in cuprates [J, H, [E S S B 0| ■ Initiated 
by the collective pinning theory 8l|, theoretical under- 
standing has been advanced [1, [jja, [13, [Hlj Ell- The 
competition between the elastic force and the random 
force builds a complex potential landscape governing vor- 
tex dynamics. An important notion has been established, 
namely the typical energy barrier felt by the system be- 
comes diver gingly large U{F) ^ at small driving 
forces 0, 0, [ill Il2j| , which results in an extremely small 
velocity of the system following the Arrhenius law. 

The degree of divergence ^ is governed by the large- 
scale elastic properties of the system and properties of 
disorder The value of exponent /i is however not 
easy to evaluate accurately, since the static and dynamic 
properties of the system are intervened in a very com- 
plex way. The situation is even more severe when the 
randomness is large and thus the order is destroyed fully. 
In most theory restoring forces assume elastic behaviors 
of the system involved, which is not the case in strongly 
disordered systems, where plastic deformations become 
important (See discussions in Ref.0]). 

We tackle this problem by computer simulations on 
three-dimensional (3D) flux lines subject to randomly 
distributed point-like pinning centers (for a parallel work 
on ID domain wall in 2D space see Refs. [Tslll^). Tuning 
the strength of random pinning force, we can reach the 
Bragg glass (BrG) [T^ . llH at weak pinning and vortex 
glass (VG) [lal at relatively strong pinning, at equilib- 
rium. Langevin dynamics then permits us to explore the 



dynamics of the system subject to driving force at various 
temperatures without any uncontrollable approximation. 

The main results are as follows: The creep law of flux 
lines predicted by previous theories is clearly reproduced. 
The exponent is estimated as = 0.5 ± 0.02 for BrG 
and /i — 0.28 ± 0.02 for VG, both universal for the re- 
spective class. While the former one is very close to the 
expected value, the latter one is a new estimation. For 
weak pinning where BrG is stable below the melting tem- 
perature Tm, we flnd an unpinned state at Tg < T < Tm 
in addition to the pinned one at T < Tg, indicating a 
replica-symmetry breaking transition at Tg ~ Tm/2. 

Model and simulation details - The model system is a 
stack of superconducting planes of thickness d with pe- 
riod s of the layer structure, with the magnetic fleld per- 
pendicular to the layers. Each plane contains vortices 
and Np quenched pins. The overdamped equation of mo- 
tion of the ith vortex at position is [13, [l8| 
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Here rj is the viscosity coefficient. The intraplane vor- 
tex repulsion is given by the modified Bessel function 
U^^ipijjZij = 0) = deoKolpij / Xab), and the interplane 
vortex attraction is U^^{pij,Zij = s) — (seo/7r)[l -|- 
ln{\ab/s)][{pij/2rgf-l] for p,j < 2rg and C/^^(py, = 
s) — (seo/7r)[l-|-ln(Aab/s)][py/rg — 2] otherwise, between 
two vortices belonging to the same flux line and sitting 
on adjacent planes, where eo = (f^ / 2t^ p^X"^^ with Xab 
the magnetic penetration depth of the superconducting 
layer, rg = js with 7 the anisotropy parameter. The 
pinning potential is U^^{pip) = —aApexp[—{pip/RpY], 
where Ap = {eod/A) ln[l + {Rl/2(i^)] with ^ab the in- 
plane coherence length and a the dimensionless pinning 
strength. Finally, F is the uniform Lorentz force, and Yth 
is the thermal noise force with zero mean and a correla- 
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FIG. 1: Velocity-force {v — F) characteristics at difFerent tem- 
peratures for a = 0.2. The solid lines are fitting curves 
based on Eq. ([2]) as described in text. The system size is 
Lj; X X A'^z = 30 X 30 X 20 with the lateral directions mea- 
sured by Xab- Errors and finite-size effects are smaller than 
the size of symbols. 



tor {FPf^{z,t)F;^f^{z',t')) = 2rjT5'P'i5{z - z')5{t - t') with 
p,q = X, y. The units of time is tq = rjX^f^/deo ~ 0.03ns 
for Xab = 2000l and ij/d ~ 10~^Pa ■ s typically for 
BSCCO [l^. Conventions including the units of phys- 
ical quantities are taken same as Ref. [isj . 

Creep motions - We start with the case of a = 0.2, 
for which an amorphous VG is realized at low tempera- 
ture [ll] . Figure 1 presents the average velocity at small 
driving forces at several typical temperatures. The driv- 
ing forces are limited to F < Fca/2 with Fco = 0.232 
the zero-temperature depinning force [isj . where the 
creep feature of vortex motion is profound. At even 
lower temperature, the dynamics is to be governed by 
the fixed poi nt at Fco, and thus omitted in Fig. 1 (see 
Refs. [Tsl. l20l|). The velocity v = 0.01 corresponds ap- 
proximately to GOm/s. 

The v — F characteristics in Fig. 1 are described well by 
the Arrhenius-type functions at respective temperatures 
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From the least-squares fitting, we find U = 0.15 ± 0.02 
and /i = 0.28 ± 0.02, with Fco = 0.232 given in the pre- 
vious work [l8| . The prefactor vq depends on tempera- 
ture roughly in a power-law way. For a stronger pinning 
a = 0.4, it is found U = 0.30 ± 0.03 and n = 0.28 ± 0.02 
(raw data not shown here). Therefore, the exponent /i is 
universal for VG. 

We then turn to the weaker pinning force a — 0.1 
where BrG is stable at low temperatures T < T™ ~ 
0.075. The v — F characteristics are displayed in Fig. 
2. Carrying on the same analysis mentioned above, we 



FIG. 2: v — F characteristics same as Fig. 1 except for a = 0.1. 



obtain U = 0.034 ± 0.006 and p ^ 0.50 ± 0.02. For 
a = 0.05, for which the equilibrium state is also a BrG 
at low temperatures, we obtain U = 0.018 ± 0.003 and 
/i = 0.50 ± 0.05. The above results indicate that the ex- 
ponent fj, for a < 0.1 falls into another class, universal 
for weak pinning strengths for which BrG is realized at 
low temperatures. 

Our estimate on the exponent fj, is in good agreement 
with the theoretical prediction 0, [ll| 
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with Ccq = the roughness exponent for BrG. 

As pointed out in Ref. 3], previous theories on vortex 
creep started from the hypothesis of elastic energy, which 
is established only for elastic deformations, and thus work 
at best for BrG. The exponent p should be reformulated 
for strong pinnings associated with VG where deforma- 
tions are plastic |3| . It was not clear yet that whether fi 
is larger in BrG or in VG [1, [2l| . The present simulation 
results suggest a stronger divergence of energy barrier in 
BrG. 

New phase boundary in BrG - So far, we have concen- 
trated on low temperatures, where the thermal energy 
only activates flux lines from pinning centers, resulting in 
an Arrhenius-type motion with a well defined energy bar- 
rier which depends on the driving force. As temperature 
is lifted, thermal fiuctuations become more important, 
which makes the competition between the randomness 
and the intervortex forces very subtle. Here we focus on 
the weak pinning case a — 0.05, for which the ground 
state is a BrG with T„i ~ 0.075. 

Figure 3 presents the v — F characteristics at several 
typical temperatures. At low temperatures, the v — F 
curves are nonlinear at low forces, similar to Figs. 1 and 
2. However, above Tg ~ 0.035, the v — F characteristics is 
linear down to the small force limit. This indicates that 
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FIG. 3: V — F characteristics for a — 0.05 at different temper- 
atures. A transition temperature Tg ~ 0.035 is defined which 
separates the nonlinear and linear v — F curves. The soHd 
lines are for eye-guide. Inset: a schematic phase diagram of 
vortex states. 



the potential barrier sensed by BrG is smeared to zero 
in an intermediate temperature regime Tg < T < T„, 
where the crystalline order is still preserved. A similar 
change is observed for relative strong pinning where VG 
is realized at low temperatures, for which the details will 
be reported elsewhere. 

We have found that even at equilibrium, i.e. for zero 
driving force, the system behaves in qualitatively differ- 
ent ways for temperatures below and above Tg. In Fig. 4, 
we display two-dimensional vortex trajectories and the 
corresponding structure factors. At T = 0.03(< Tg) 
(Fig. 4(a)), vortices are trapped in cages formed by the 
random pinning potential and the intervortex repulsions 
(even for time much longer than that shown in Fig. 4); a 
profound BrG order is clearly seen. At T = 0.09(> T„i) 
(Fig. 4(c)), the system is in a liquid state where vortices 
diffuse freely and randomly. At T = 0.06, the behaviors 
of the vortices are different as evidenced in Fig. 4(b): 
Over an intermediate time scale, vortices are trapped at 
local positions in a way such that the BrG order is es- 
tablished. Vortices then move quickly to another set of 
localized positions, which also establishes the BrG order, 
and stay in the new positions over another intermedi- 
ate time duration. The intermittent trapped-and-moving 
motions continue, and eventually vortices can travel over 
distances larger than the vortex lattice constant qq . 

In order to capture the vortex motions better, we have 
monitored the number N{t) of vortices that move over 
distances smaller than ao after time t. As shown in 
Fig. 4(d), N{t) decreases exponentially with t in a large 
time scale at T = 0.06 since vortices move over distances 
larger than the lattice constant. The exponential decay 
occurs even in a small time scale at T = 0.09 as the sys- 
tem is random. At T = 0.03, N{t) remains unity during 
the time evolution because all the vortices are caged. 

We have also studied how vortices change their near- 




FIG. 4: Trajectories of vortices on a typical layer and struc- 
ture factors averaged on all layers at three typical tempera- 
tures: (a): T = 0.03, (b): T = 0.06, and (c): T = 0.09. Only 
a part {Lx x Ly — 14 x 14) of the totally simulated system is 
displayed. The time interval is lOOro and the total time shown 
here is 50,000ro. (d): temporal evolution of the number of 
vortices which are trapped in the initial cages. 



est neighbors (NNs). At t = 0, we identify the NNs 
of each vortex in terms of the Delaunay triangulation 
method. Due to random motions, some vortices defined 
as NNs a.t t — diffuse away. The quantity CNN(i) is 
defined as the number of NNs identified at f = remain- 
ing as NNs up to time t, with a trivial normalization 
to unity at t — 0. The quantity Cnn varies with time 
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FIG. 5: Temporal evolution of number of nearest neighbors 
Cnn kept from the initial state at t = (see the definition in 
text) for Q = 0.05 at three typical temperatures. 
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FIG. 6: Fluctuations of the structure factors at the first 
reciprocal lattice vector S{Q) (a) and fraction of sixfold- 
coordinated vortices Pe (b) for a = 0.05 at three typical 
temperatures. Green, blue and red are for T = 0.03, 0.06, 
and 0.09, respectively. 



in qualitatively different ways in the three temperature 
regimes as presented in Fig. 5: At T = 0.03, Cnn re- 
mains unity as time evolves where vortices are trapped 
in the initial cages and thus NN's remain the same. At 
T = 0.09, Cnn decreases to zero exponentially fast be- 
cause the system is totally random. At T = 0.06, Cnn 
decreases slowly with time. Checking the trajectories of 
vortices, we find that, in the intermediate temperature 
regime Tg < T < Tm, diffusions of vortices take place in 
terms of swapping nearest neighbors. Due to the mobil- 
ity of vortices at Tg < T < Tm at equilibrium, the v — F 
characteristics should be linear even at infinitesimal driv- 
ing force as revealed in Fig. 3. 

Along with the random motions and neighbor swaps of 
vortices, we also monitor the variation of crystalline or- 
der of the whole vortex system. In Fig. 6, the structure 
factor S{Q) at the first reciprocal lattice vector of the 
triangular lattice and the fraction of sixfold-coordinated 
vortices Pg are displayed for the same time duration of 
Fig. 5. For T = 0.06, the value of S{Q) remains con- 
stantly ~ 0.3 and Bragg peaks can be clearly observed 
(see Fig. 4(b)). Therefore, the system always keeps the 
ordered triangular structure with most vortices sixfoldly 
coordinated, even though the individual vortices are mo- 



bile in a stochastic way. Compared with thermal fluctua- 
tions, vortex- vortex interactions are strong while pinning 
is weak, resulting in an unpinned BrG. The pinned BrG 
at T = 0.03 is well ordered as captured by large values 
of S{Q) and Pq, while they are both small at T = 0.09 
corresponding to the vortex liquid. 

Recently an experimental finding of a second-order 
phase boundary within the solid phase of the vortex sys- 
tem in BSCCO was reported [22|; the transition line 
Bg{T) dividing the BrG phase domain into two parts. 
This transition was discussed as a transition between 
the two different types of glasses, the high-temperature 
pre-glass or marginal glass phase, with the nearly linear 
response behavior, and the low-temperature true glass 
domain. The transition between the two phases man- 
ifests itself in a replica-symmetry breaking within the 
replica description of disordered vortex system [23]. It 
is tempting to associate our observation of the change 
in dynamic vortex behavior with this transition. Indeed, 
at low temperatures T < Tg, vortices remain indefinitely 
near their equilibrium positions since on top of the elas- 
tic forces keeping them there, vortices are additionally 
immobilized by disorder. At intermediate temperatures, 
Tg < T < Tm, pinning of single vortices is not efficient 
any more and they may switch between their equilibrium 
positions due to thermal diffusion. However, putting this 
conclusion on a firm quantitative basis requires a more 
detailed study of the processes of vortex relaxation in 
both phases and will be a subject of a forthcoming pub- 
lication. 
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